function [ w ] = wendlandWeight( d, rho )
%WENDLANDWEIGHT Compute the wendland weight function with support [0, rho]
%   d   : distance between the xi node and the sample point x
%   rho : upper bound of the function support

w = (1-d/rho).^4.*(4*d/rho+1);

end

